How hard is it to approximate the Jones polynomial? 
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Freedman, Kitaev, and Wang floll . and later Aharonov, Jones, and Landau 0], established a quantum algo- 
rithm to "additively" approximate the Jones polynomial V(L,t) at any principal root of unity t. The strength 
of this additive approximation depends exponentially on the bridge number of the link presentation. Freedman, 
Larsen, and Wang fl^ established that the approximation is universal for quantum computation at a non-lattice, 
principal root of unity; and Aharonov and Arad 0] established a uniform version of this result. 

In this article, we show that any value-dependent approximation of the Jones polynomial at these non-lattice 
roots of unity is #P-hard. If given the power to decide whether |K(L,r)| > a or |K(L,f)| < b for fixed constants 
a> b>Q, there is a polynomial-time algorithm to exactly count the solutions to arbitrary combinatorial equa- 
tions. In our argument, the result follows fairly directly from the universality result and Aaronson's theorem 
that PostBQP = PR d. 



1. INTRODUCTION 

A well-known paper of Aharonov, Jones, and Landau ^ 
establishes a polynomial quantum algorithm to approximate 
the Jones polynomial at any principal root of unity; a more 
abstract form of this algorithm appeared previously in a paper 
of Freedman, Kitaev, and Wang lUOll . 

Theorem 1.1 (Freedman, Kitaev, Wang [10]; Aharonov, 
Jones, Landau Q)- Let t = exp(2;r//r), let L be a link pre- 
sented by a plat diagram with bridge number b, and letV{L, t) 
be its Jones polynomial. Then there is a polynomial-time 
quantum decision algorithm that answers yes with probability 



P[yes] 



V{L,t) 



(fi/2 + f-i/2) 



b-l 



(In the AJL version of the result, the algorithm is jointly 
polynomial time in the r, the order of the root of unity; as well 
as in the bridge number and the crossing number They also 
refine the algorithm to estimate V{L,t) as a complex num- 
ber rather than just estimating its length.) Aharonov, Jones, 
and Landau describe the error in this algorithm as additive, 
and note that it would be much harder to provide algorithm 
with multiplicative error Multiplicative error would mean that 
V{L,t) or |y(L,f)| can be approximated to within a constant 
factor of its value. 

Another way to distinguish between types of error is to say 
that the error in the approximation algorithm is presentation- 
dependent. Given different plat presentations of the same link 
L, the error grows exponentially in one of the parameters of 
the presentation, namely the bridge number. An algorithm to 
approximate the Jones polynomial is only directly useful for 
topology if its error is value-dependent, i.e., independent of 



quantities other than the value of | V(L, f ) |. Multiplicative error 
is one type of value-dependent error, but in some cases value- 
dependent error is more general. (Multiplicative error is used 
to define complexity classes such as APX; see Section lTTl ) 
Freedman, Larsen, and Wang lfl2ll established that the 



approximated quantity 



in Theorem 11.11 is 
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universal for quantum computation when r = 5 or r > 7. 
(Aharonov and Arad 10] establish an r-uniform version of this 
result.) Thus, even if the error is presentation-dependent, the 
approximation is computationally valuable for carefully cho- 
sen link diagrams. 

On the discouraging side, Jaeger, Vertigan, and Welsh ITtIi 
showed that it is #P-hard to exactly compute the Jones poly- 
nomial V{L,t) except when f is a root of unity of order 
rG {1,2,3,4,6}. They established a reduction from the Tutte 
polynomial of a planar graph to the Jones polynomial of an 
associated link. 

The main result of this article is that the "encouraging" 
universality result strengthens the "discouraging" hardness re- 
sult: Any value-dependent approximation of above values of 
the Jones polynomial is #P-hard. The argument is a mash- 
up of three standard theorems in quantum computation: The 
Solovay-Kitaev theorem ll20ll (although only a mock version is 
needed), the FLW density theorem, and Aaronson's theorem 
that PostBQP = PP H. (See also H for a different hardness 
result.) 

Theorem 1.2. Let I{L) be a quantumly finite, unitary link in- 
variant which is eventually dense. In particular, I{L) can be 
the Jones polynomial V{L,t) with t = exp(2;r//r) and r = 5 
or r >1. Let a > b > be two fixed positive real numbers, 
and assume as a promise that \I{L)\ > a or \I{L)\ < b. Then 
it is #P-hard, in the complexity of the diagram ofL, to decide 
which inequality holds; moreover we can take L to be a knot. 

The idea of Theorem 1 1 .21 applies to many cases where the 
invariant is not unitary; and when it is not an invariant of links 
but some other type of topological or combinatorial object. 
We have no formal statement of a general result, but the ba- 
sic argument is that if an invariant can simulate a quantum 
computer, then typically multiplicative or value-dependent ap- 
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proximation is universal for PostBQP and therefore #P-hard. 
Here is an example result of this type, modulo a denseness 
result which is explained in Section |4] and which will be es- 
tablished in a later paper (We also provide a Sage program to 
confirm an adequate weak form of the conjecture for specific 
values of n.) 

Theorem 1.3. Multiplicative approximation to the Tutte poly- 
nomial T{G,x,y)for planar graphs G is #P-hardfor any fixed 
values ofx andy such that n — (x~l){y—l) >4 andx,y < 0, 
except X ~y ^ —1. Multiplicative approximation of the Jones 
polynomial is also #P -hard for any real f > other than t = \. 

In related results, Goldberg and Jerrum Iil4l1 showed that 
multiplicative approximation of many values of the Tutte 
polynomial T{G,x,y) is NP-hard (relative to RP) for non- 
planar graphs, while some values are #P-hard. Jaeger, Verti- 
gan, and Welsh UtIi also analyzed when T{G,x,y) is #f-hard 
to compute exactly. They noted that the Jones polynomial 
y(L,f) of an alternating link L is equivalent to T{G,x,y'\ for 
a planar graph G along the curve xy = 1. More recently lUSll . 
Goldberg and Jerrum also established that many values of the 
planar Tutte polynomial are NP-hard to approximate. Their 
new theorems have some intersection with Theorem II. 31 but 
their constructions are very different. 
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2. COMPLEXITY THEORY 

We assume that the reader is somewhat familiar with com- 
plexity classes such as P, NP, BQP, #P, and the notation that 
means the class A with oracle B. See the Complexity Zoo 
llzsll and Nielsen and Chuang L20i1 for a review. 

Whereas a problem in the class #P counts the number of 
witnesses accepted by a verifier in polynomial time, and a 
problem the class N P reports whether there is an accepted wit- 
ness, a problem in the class PP reports whether a majority of 
the witnesses are accepted. A problem is #P-hard if and only 
if it is PP-hard, because P'''^ = P*'' (exercise for the reader). 
On the other hand, a problem which is #P-hard is also hard 
for the polynomial hierarchy PH, by the deeper theorem due 
to Toda lUl] that 

oo NP 

ph'=i'Unp^^cp#p. 

"=1 n 

The class NP with a tower of « — 1 NPs as an oracle is called 
the «th level of the polynomial hierarchy. One of the standard 
conjectures in complexity theory is the polynomial hierarchy 
does not collapse, i.e., that nth level does not equal the « + 



1st level for any n. Thus by Toda's theorem, if a problem is 
#P-hard, then it is viewed as qualitatively harder than if it is 
merely NP-hard. 

The class PostBQP is defined as polynomial-time quantum 
computation with free retries, or postselection. In other words, 
the computation can output "yes", "no", or "try again", and 
the probability of "yes" is required to be c times that of "no" or 
vice-versa for some constant c > L All values of c are equiva- 
lent because c can be amplified by repeated trials. There is an 
analogous class PostBPP for classical randomized computa- 
tions; it was also defined previously as BPPpath- Aaronson [Ql] 
estabHshed that PostBQP = PP. It is not hard to show that 
PostBQP is a subset of PP, just as BQP, NP, and a number 
of other important classes are. The more surprising fact is that 
PostBQP is all of PP. 

By contrast, PostBPP is unlikely to be all of PP. The rele- 
vant complexity results are as follows: 

1. PostBPP contains pH^iP (P with parallel NP queries) 113]. 

2. PII'^P equals P^pM (P with logarithmically many NP 

queries) 

3. PostBPP derandomizes to pH'^^, i.e., they are equal if suf- 

ficiently good pseudo-random number generators exist 

iH. 

4. Without any derandomization assumption, PostBPP C 

BppNPcNP^P"" 13. 

Thus, PostBPP is known to be in the third level of PH and it 
is presumably in the second level. 

Another interpretation of PostBQP or PostBPP is that it 
consists of two computations performed by Alice and Bob. 
Given a decision problem L, Alice and Bob each either suc- 
ceed or fail in separate polynomial time with probabilities a 
and b. If L{x) is "yes", then a > cb, while if L{x) is "no", 
then b > ca. This description already reveals the main idea of 
Theorem lL2l If we can estimate both a and b with error less 
than c, then we can compare them up to a factor of c^. 

It is interesting and relevant to compare PostBQP and 
PostBPP to three other complexity classes: AgPP (one-sided 
almost wide PP), SBP (small-bounded probabilistic P |5i|), 
and what we will call SBQP. All three classes depend on 
a real-valued function f{x) in FP (expressed in fixed-point 
arithmetic, say), where x is the input to the decision problem, 
and a constant c > I. The classes SBP and SBQP are defined 
in exactly the same way as PostBPP and PostBQP, except 
that Bob's value b is not a success probability but rather com- 
puted directly in FP (as a real number in fixed-point arith- 
metic, say). Meanwhile AqPP is defined like SBP, in which a 
is a success probability and b is computed in FP, except with 
the formulas 

L{x) — yes =^ a> cb + — 

1 1 

L(x) ~ no ==> - < fl < D + - . 

Proposition 2.1. Without loss of generality, the threshold 
function f(x) in the definition o/ AqPP, SBP, SBQP can be 
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taken to be 2 ''^'^l^ for some polynomial p; and all values of 
the constant c are equivalent. 

Proof. The constant c is irrelevant by the usual technique of 
amplification by repeated trials. This is immediate in the case 
of SBP and SBQP. It is not very difficult in the case of AqPP, 
and was established by Vyalyi Ii24l1 . 

To argue that f{x) can be set to 2^''(l-^l) (in the cases of 
SBP and SBQP), first choose p so that f{x) > 2-''^\^\\ Then 
Alice can compute f{x) and reduce her success probability by 
a factor of 2''(l'''l^/(x). The argument in the case of AqPP is 
essentially the same and was also explained by Vyalyi ll24ll . 



If we assume that b > j and let C = 2 be the threshold factor 



□ 



Proposition 2.2. 



SBQP = AoPP. 



Proof The proof is almost the same as Aaronson's proof that 
PostBQP = PP. We can also define AqPP as a counting class 
in which, for each certificate y of length «, the computation 
produces a value /(y) = ±1, and these values are summed to 
produce A (;ic). For a decision problem L e AqPP, we require 
that 



L{x) ~ yes 
L{x) = no 



A{x)>2"Cb 
0<A{x)<2"b. 



First, let L G SBQP be computed by a quantum circuit that 
consists of Hadamard and Toffoli gates. It is convenient to 
change the counting model of AqPP slightly to let the values 
be ±1 or 0. Then we obtain an AqPP algorithm by multilin- 
ear expansion of the effect of these gates on density matrices. 
The matrix entries of a Toffoli gate, in its effect on a den- 
sity matrix, are and 1; the corresponding matrix entries of 
a Hadamard gate are ±5. The final probability is given by a 
partial trace of the output density matrix, and is non-negative 
and exactly matches the criteria for AqPP. 

Now let L e AqPP and let a be Alice's success probability 
in the AqPP algorithm. We can again slightly re-express the 
counting model of AqPP so that f{y) G {0,1} and its sum 
A ~ A{x) is given by A = 2"fl. 

Then, in the SBQP algorithm, we can quantum-compute 
the unitary map 

where the value f(y) is written to an ancilla qubit. We provide 

the input | + H h) to JJf, and then postselect on whether 

the left n qubits of the result are all |+). If they are, then the 
ancilla qubit has the state 

|V/)-(l-fl)|0)+«|l). 

If this qubit is measured in the ± basis, then the probability of 

|-)is 

, (2fl-l)2 



in the AqPP algorithm, then 



1 



< fl < - +/? 
2 



a > - + 2/7 



a < 



a > 



I+4Z72 
16^2 



1 + 



So we can let b' ~ 4-b^ and C' = 2 in an SBQP algorithm that 
produces the probability a'. □ 

Many of the complexity classes discussed here employ the 
semantic condition that the probabilities of particular out- 
comes are above or below certain thresholds. We can also con- 
sider promise versions of these classes in which this threshold 
gap may only present for some inputs and not others. When 
they are considered in promise form, SBQP and SBP become 
equivalent for our purposes to PostBQP and PostBPP. 



Proposition 2.3. 



l + (2fl-l)2- 



PostBPP C pP'-o'Tii^^SBP 
PostBQP C pPromiseSBQP 



Proof Suppose that L G PostBQP (say) and that it is imple- 
mented by a quantum circuit. We assume that the gates of 
the circuit are Hadamard and Toffoli gates (see the remark be- 
low), so that the probabilities a and b, if they are non-zero, 
are bounded below by 2^'^'*'' in a circuit with g gates. They 
are therefore also bounded below by 2^P°'y'l*l', where x is the 
input. 

The construction is then similar to another part of Aaron- 
son's proof that PostBQP = PP. We assume that either a>%b 
or that b > 8a. Choose an n, no more than polynomial in 
such that max{a,b) > 2^". Then for each < A: < n, use 
PromiseSBQP to compare both a and b to 2^*^. If a > 8fe, 
then for every k, PromiseSBQP will either reliably report that 
a > 2^* or that 2^^ > b, and there will exist a k for which it 
will do both. Meanwhile if b > 8fl, it will report that b > 2^* 
or that 2^'^ > a, and both for at least one k. These two out- 
comes are mutually exclusive. 

The argument that PostBPP C pPromiseSBP ^^^^^ 
simpler since the lower bound on ma\{a,b) is immediate. □ 

Remark. Gate-set independence is a non-trivial property of 
PostBQP. Because the probabilities a and b are polyno- 
mials in the matrix entries of the gates, and using norm 
bounds on the coefficients of those polynomials, the gates 
can be varied slightly so that a,b > 2^?°'^*^^'. The Solovay- 
Kitaev theorem can then be applied to show that PostBQP 
is a gate-independent complexity class. By contrast, we do 
not know whether post-selected quantum computation is gate- 
independent with a time bound of (5(«") for some fixed expo- 
nent a, because the Solovay-Kitaev theorem could change the 
exponent. 

Finally, as noted by Aaronson, linear computation is 
another interesting interpretation of PostBQP. Post- 
conditioning allows us to replace unitary gates by subunitary 
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gates, and to rescale subunitary gates arbitrarily. But every 
linear operator that acts on vector states | is proportional to 
a subunitary operator. Thus, PostBQP can also be defined by 
the class of polynomial-sized circuits with linear gates, with- 
out the unitary restriction. It may seem then as if PostBQP 
uses the state space of a qubit only as a 2-dimensional com- 
plex vector space, and not as a Hilbert space. However, the 
definition of measurement probabilities in PostBQP still use 
the Hilbert space structure. 



2.1. Approximation classes 

The approximation classes listed in the Complexity Zoo 
llzsll that express multiplicative error include APX, PTAS, 
and FPTAS. These classes are defined there for optimiza- 
tion problems, but they can equally well be defined for arbi- 
trary functional problems. Let f{x) be a function f{x) that 
takes strings x to positive real numbers. Then f{x) is in APX 
if it can be approximated to within some bounded factor in 
polynomial time (with fixed-point output); it is in PTAS if it 
can be approximated to within a factor 1 + e in polynomial 
time for any e > 0; and it is in FPTAS if the computation is 
jointly polynomial time in |x| and 1 /e. (These classes all refer 
to deterministic computation; there are analogous randomized 
classes such as FPRAS.) 

We do not know of a standard complexity class to express 
general value-dependent error, so we define one here. Let f{x) 
be as before. Then / is in APV if for every constant a > 0, 
there exists a constant b> a and a polynomial-time algorithm 
to decide whether f{x) > boi f{x) <a, given the promise that 
one of the two is true. Similarly, we could define a randomized 
version ARV; and both APV and ARV have a variation which 
is uniform in the choice of the constant a. 

The following proposition says that if f{x) can be suitably 
rescaled, then general value-dependent error becomes equiv- 
alent to multiplicative error. The proposition and its proof 
use the same idea as the proof of Proposition |23J and a sim- 
ilar rescaling argument in Aaronson's argument that PP C 
PostBQP. 

Proposition 2.4. Suppose that f{x) takes positive real values 
and is in APV, and suppose further that \ log/(jc) | is bounded 
by a polynomial in the length \x\ of the input. Suppose that 
there are constants c > 1 and k> \ such that every integer n, 
there is a reduction yn{x) such that 

f{x)<k"f{y„(x))<cf{x), 

and suppose that this reduction can be computed in joint poly- 
nomial time in n and in \x\. Then f{x) is in APX. 

Proof. (Sketch) Let a and b be some constants such that we 
can decide by a subroutine whether f{x) < a or f{x) > b in 
polynomial time. Then we can bound f{x) to within a factor 
of cb/a. We know by hypothesis that f{x) > k^'" and f{x) < 
k"^ for some m which is polynomial in \x\. So the strategy is 
to ask whether f{yn{x)) is more than a or less than b for every 
«| < m. The largest « for which the subroutine reports that 



f{yn{x)) < a yields a good estimate of ak^". The estimate is 
within a factor of cb/a, even though the subroutine could give 
a false yes answer when f{yn{x)) < b. □ 

3. QUANTUM INVARIANTS 

Let A: be a field and let I{L) be a A:-valued invariant of links. 
Assume that I{L) is multiplicative, meaning that if L = Li UL2 
consists of two links separated by a plane, then 

/(L)=/(Li)U/(L2). 

Depending on the nature of /, the link L could optionally have 
other decorations such as a framing, a coloring, or an orienta- 
tion. (An orientation refers to orienting each component of the 
link as a curve; ambient is always oriented by convention.) 
Now choose a fixed plane C and consider those links 
L ~ AUB that are formed from two tangles A and B. Then 
/ can also be seen as an invariant I{A,B) of the two tangles, 
defined when A and B have matching boundary E so that they 
can be glued together to make a link: 




We thus also have a relative invariant /(A) whose 
value is a function of the missing tangle B. Let Tz be the set of 
all possible choices of B. Then /(A) technically takes values in 
the infinite-dimensional vector space of functions Fun{£^,k). 
But it may happen that the span of the values of I (A) is 
a finite-dimensional subspace of Fun{3^,k). The invariant / 
is called quantumly finite if it has this property for all choices 
of E. Note that / induces a non-degenerate bihnear form on 
y(E). Indeed, / is a bilinear form on Fun(5^,^) and V(E) is 
the quotient of Fun (5^,^) by the kernel of /. 

(For related constructions, see also 10] and Jgt].) 

A linear relation among tangles in J^, whether or not it is 
finite-dimensional, is called a skein relation. The space 
itself is called a state space or a skein space. Also, the mul- 
tiplicative property or L implies that if E is divided into two 
disjoint parts, E = Ei UE2, then 

y(Ei)(8)y(E2)'^=^. 

The state space is also naturally a representation of the 
braid group of E. In the case of the Kauffman bracket, this 
can be taken to be the usual braid group B„, but in general it 
will be a braid group that preserves the colors or orientations 
of the points in E. 
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If A and B are two matching tangles as before, let A and B 
be their reflections through their separating plane. If the field 
^ = C and / conjugates under reflection, then 



I{B,A)=I{A,B). 

In this case I{A,B) induces a non-degenerate Hermitian form 
on If this form is positive-definite, then is naturally 
a Hilbert space and the invariant / is unitary. If / is unitary, 
then we can also write it and the associated tangle invariant in 
Dirac notation: 

{L)j=I{L) \A)j=I{A). 

Say that a unitary link invariant / is fully non-trivial if it 
satisfies the following two conditions. 

1. For every pair of colors, left and right crossings are linearly 
independent in the their skein space 



2. If Uc is an unknot with any color c, then {uc)i > 1. 

In fact, the second condition follows from the first one for 
unitary invariants, but it is simpler for us to take both con- 
ditions as the definition. Note in particular that if / is fully 
non-trivial, then A\m,y^ > 2. 

An important example of these axioms is the Kauffman 
bracket form of the Jones polynomial. The bracket is defined 
for unoriented links with "blackboard framing", meaning that 
the link is a ribbon laid flat along the link diagram. If f is the 
indeterminate of the Jones polynomial, then the defining skein 
relations are 



= -fV4| 



-1/4 



K 
,-1/2 



(Note also that the Jones polynomial is often also parameter- 
ized by the variable q^t^ .) Let V{n) be the Kauffman state 
space of tangles with n endpoints. Then the Kauffman rela- 
tions imply that 



dimy(2«) = C„ = 



1 



1 



the nth Catalan number, except when f is a root of unity, 
whence it is can be smaller On the other hand, if f = e^"'/'' is 
a principal root of unity of order r, then the Kauffman bracket 
is a unitary invariant of framed links. The Jones polynomial 
is defined from the Kauffman bracket in turn by the relation 

jmt)^{~t'r<^\L)K. 

Here w{L) is the writhe of the diagram of L, i.e., the number 
of positive crossings minus the number of negative crossings. 
(The Jones polynomial is also quantumly finite after rescaling 
to make it multiplicative, but the Kauffman bracket formula- 
tion is generally more convenient.) 



If I{L) is a quantumly finite invariant of links, then there is 
an associated action of a braid group on each skein space j/^. 
If E has « equivalent marks, as in the case of the Kauffman 
bracket, then we can take the braid group to be the standard 
braid group B„. (In more general cases, the braid group has to 
preserve the colors or orientations of the marks on E.) If I{L) 
is unitary, then we say that it is eventually dense if the image 
of the braid group is dense in PSU(c;5^) when E has enough 
points (of every available color, if there is more than one). 

Theorem 3.1 (Freedman, Larsen, Wang lfl2ll '). The Kauffman 
bracket is eventually dense ift = e^^'^'' with r = 5 or r >7. 

The following proposition was first argued by Freedman, 
Larsen, and Wang in a special case tfTT tl . Our argument is not 
really different from the standard one, but we give it for com- 
pleteness. For simplicity, we assume that / is an uncolored 
link invariant in the rest of the section. 

Proposition 3.2. Suppose that a quantumly finite, fully non- 
trivial link invariant I is eventually dense and that for each 
fixed L, the value of (L); has an FPTAS. Let p{x) be the prob- 
ability that some quantum-polynomial time algorithm accepts 
an input X. Then the input x can be encoded as a link L = Lix) 
with bridge number g, so that 



p{x) 



where "k," is in the FPTAS sense and u is an unknot. 



(1) 



The condition that (L)/ has a FPTAS for each L is fairly 
mild. It says that (L)/ can be approximated to within e in 
polynomial time in 1 /e, and that the polynomial can get worse 
as L gets bigger. In the proof this estimation property will 
be combined with the Solovay-Kitaev theorem. However, 
the Solovay-Kitaev theorem is much stronger than necessary: 
it establishes that quantum gates can be approximated with 
words of size polynomial in log e, and that the words can be 
found in the same time complexity. 



Proof. There is more than one way to argue Proposition 13. 21 
in this article we will use the method of plat presentation. For 
simplicity we will describe the case of an uncolored, unori- 
ented link invariant; the more general case works the same 
way. 

A plat presentation of a link L is a diagram that consists of 
a braid B on 2g strands, with n copies of the cap v at one end 
and g copies of the cup v at the other end: 
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The number g is the bridge number of the plat presentation. 

First, using uncomputation, we can make a quantum circuit 
C = C{x) so that the acceptance probabiUty is 

l(0"|C|0")|2, 

and such that both the gate number |C| and the qubit number 
n are polynomials in 

Second, by hypothesis the space V{4-) for the invariant / is 
at least 2-dimensional. We can use a 2-dimensional subspace 
of y(4) as a computational qubit. More precisely, we choose 
this subspace to include the value of 2v, where v is a cap: 



The computational state |0) is defined to be proportional to 
|2v) in the computation. Since / is eventually dense, there ex- 
ists an n such that the closure of the braid action on V{4n) in- 
cludes 2-qubit gates acting on V(8). Thus the circuit C can be 
approximated by a braid B on 4« strands. The FPTAS hypoth- 
esis implies that the matrix entries of the braid generators can 
be approximated in polynomial time, and the Solovay-Kitaev 
theorem implies that B can be generated as an FPTAS of C. 

Let L be the link which is the plat closure of B. Then 



(0"|C|0") 



where as usual u is an unknot. We let g^ln. Since this norm 
of the left side is the acceptance probability, this completes 
the proof of the proposition. □ 

In light of the analysis in Section|2]of PostBQP, Proposi- 
tion l3.2l establishes that if I{L) is unitary and eventually dense, 
then a multiplicative approximation to it is PostBQP-hard. By 
Aaronson's theorem, it is therefore also #P-hard. 

One remaining task to prove Theorem 11.21 is to refine the 
construction to general value-dependent approximation rather 
than multiplicative approximation. Explicitly, let a > ^ > be 
constants as in the statement of Theorem ll.2l let p and c be the 
polynomial the constant in the modified definition of SBQP 
in Proposition 12. II By that proposition and equation ([TJ, it is 
SBQP-complete and therefore #P-hard to determine whether 

mi? 



is more than c2 ^'1*1' or less than 2 ^'1' 
#P-hard to determine whether 



In other words, it is 



mi) 



2 J >c|(m)/|«2-''(W) 



< c2-P(l- 



We want to make a modified link L' to make it hard to deter- 
mine whether | (L') ^ is more than a or less than b. Recall that 
i> is a polynomial in |a| and that | («)/| > 1. If 

|(m)/|'^'«2P 



when |.y[ is large, then we can add m{\x\) copies of the unknot 
to L so that 

\{u)i\<^+"'2-P 
is bounded. On the other hand, if 

2(log2|(«)/|)«(W)»p(|x|) 

then we can use denseness to first create a link Lq (say a 2- 
bridge link corresponding to a 1 -qubit circuit) such that | (Lo) | 
is a small constant. Then we can add ?m(|x|) copies of Lq to L 
so that 

i(io)n(«)/i*'2-'^ 

is bounded. The constant c in the definition of SBQP can be 
chosen to overwhelm the bound in either case as well as the 
specific values of a and b. 

Finally we want to show that L can be a knot. The trick for 
this is that since the braid group is eventually dense, the pure 
braid group is also eventually dense. Thus in the braid B we 
can switch two strands, and then approximately cancel its ef- 
fect with a pure braid that does not permute any strands. The 
permutation induced by B is thus decoupled from the approx- 
imate value of (L)/, so it can be chosen so that L has only one 
component. This completes the proof of Theorem |1.2l 

4. THE TUTTE POLYNOMIAL 

The aim of this section is to prove Theorem I L3I In order 
to understand the Tutte polynomial, we will first convert it to 
another graph invariant with equivalent information known as 
the Potts model. The Potts model of a graph G depends on a 
positive integer n, the number of colors; and on a variable y. 
The weight of a coloring of the vertices of G with n colors is 
defined as y*^ if k of the edges of G connect two vertices of 
the same color Then the Potts partition function Z{G;n^y) is 
defined as the total weight of all vertex colorings. The Potts 
partition function yields the Tutte polynomial T{G\x,y) by the 
formula 



T{G;x,y) = (y - - l)-'WZ(G;«,y), 

where n ^ {x — \ ){y — \ ). A more standard definition of the 
Tutte polynomial is via a contraction-deletion formula. How- 
ever, the definition using the Potts model is easily equivalent. 

In order for the complexity of a value of the Tutte polyno- 
mial or Potts model to be well-defined, the parameters them- 
selves must be well-defined computationally. One conserva- 
tive choice is to let x, y, and n be rational numbers. However, 
for our purposes we can let x, y, and « be any real numbers that 
are themselves defined by an FPTAS. (Or complex numbers, 
although this section has no bearing on them.) 

We can define a graph with boundary to be any graph G to- 
gether with a marked subset dG of the vertices, by definition 
the boundary vertices. Given two graphs A and B with bound- 
ary and a bijection between dA and dB, we can glue them 
together to make a graph G, in the same way that two tangles 
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can be glued together to make a link. Thus for any invariant 
I{G) of graphs, we can also define the relative invariant I{A) 
and say what it means for / to be quantumly finite. There is a 
state space V (k) for each number of vertices k. We can also 
consider the same constructions for planar graphs; one differ- 
ence is that the boundary vertices dA of a planar graph A must 
all lie on the outside face of A. 

Although we could define the Potts model Z and the state 
space V{k) over any field, in this section we will work over 
the real numbers R. 

Proposition 4.1. Ifn and y are indeterminates, then for the 
Potts model for all graphs, 

A\mV{k)^Bk, 

where Bi^ is a Bell number, the number of partitions of a set 
with k elements. For the Potts model for planar graphs, 

dimy(fc) 

the kth Catalan number 

Here is an example of a shrub forest, in this case a planar 
one, which can be used as a basis vector for V{k): 




I 



Proof (Sketch) This is a folklore result in enumerative com- 
binatorics and quantum algebra. For each partition P of the 
boundary k vertices, we can construct a shrub forest such that 
each parts of P is the set of leaves of one of the shrubs. In the 
planar case, we can consider planar shrub forests and planar 
partitions. In both cases, the shrub forests are a basis ofV{K). 

The fact that shrub forests, respectively planar shrub 
forests, span follows from the contraction-deletion formula 
for the Potts model. The fact that they are linearly indepen- 
dent can be established by writing down the matrix I{A,B), 
where A and B are two shrub forests. For each entry of this 
matrix, look at the leading term as a polynomial in y, and then 
look at the leading term in « of its coefficient. It is not hard to 
show that this leading term of the diagonal term of det/(A,B) 
cannot be cancelled, so that I{A,B) is non-singular. □ 

An important variation of the Potts model (or the Tutte 
polynomial) is the multivariate version, where the weight y 
can be different for each edge of G. It will be important to as- 
sign each edge not only its weight y, but also the dual weight 
X using the relation n = (y — 1 ) (x — 1 ). (Note that there is only 
one value of n for all of G.) 

The Potts model admits a renormalization procedure so that 
the ordinary Potts model becomes equivalent, for the purpose 



of #P-hardness of multiplicative approximation, to part of the 
multivariate Potts model liT4ll . If G has two parallel edges 
with weight yi and y2, then they are equivalent to a single 
edge with weight yiy2- Meanwhile, if G has two edges in 
series with dual weight xi and X2, they are equivalent (up 
to changing Z by a constant factor) to one edge with weight 
xiX2. These transformations are called shift operations; they 
are also called compositions and implemented weights lfl3ll . 
Note that the series and parallel compositions also preserve 
planarity. 

The point of implemented weights is the following reduc- 
tion: 



Proposition 4.2. Let POTTS (Fjn) be the problem of comput- 
ing the Potts model with n colors and with a set Y of allowed 
weights. Let Y' be all weights that can be implemented from Y : 
the closure ofY with respect to multiplication of the weights y, 
with respect to multiplication of the dual weights x, and with 
respect to limits of sequences of weights with polynomial over- 
head. Then multiplicative approximation of ^OYY^iY' ^n) re- 
duces to multiplicative approximation o/POTTS(y,n). 



Goldberg and Jerrum ifTill establish the following case of 
Proposition l4.2l 



Proposition 4.3. Consider graphs with the Potts model with 
n colors and with y {i.e., Y = \y\). Suppose that n > 2 and 
x,y <0, except for the pointy = x = — I. Then all weights y' ^ 
1 can be approximated with polynomial overhead by parallel 
and series compositions. 



Proof (Sketch) First, for any n > 1, suppose that we can ob- 
tain some weight y\ and let x' be the dual weight. Then by 
taking powers of the dual weight x', we can obtain weights 
y that converge to 1 . By taking powers of these weights, we 
can obtain a dense set of weights /" > 1 . With more careful 
accounting, this strategy converges to any weight y'" > 1 with 
polynomial overhead. 

Now suppose that « > 2 and that y>—l. Then y' =y^ > 1 , 
so by the previous argument, we can obtain all weights y" > 1 . 
We can now multiply these weights by the original weight y, 
and the dual weights x" by the dual weight x, to obtain all 
weights >'"' < 1. Likewise if jc > ~1, we can also obtain all 
weights. 

Suppose instead that x,y > — 1. By hypothesis, either 
0>x> — 1 orO>3'> — 1; suppose the latter without loss 
of generality. Then we can let / = y" for a large enough odd 
power n, and obtain x' < — 1 . Thus we can reduce to the pre- 
vious case. 

The following diagram shows the different cases for two 
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values of n. 
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We can define two types of operators that act on the state 
space V{k) and that generate a semigroup which is analogous 
to the braid group action on the state space of the Jones poly- 
nomial in Section [3] The operators are expressed by graphs 
with two subsets of vertices, which may intersect, and which 
are called left and right boundary. One type of operator Ay 
is given by connecting two of the vertices which are adjacent 
with an edge with weight y; all k of the vertices are both left 
and right vertices. The other type of operator is given by 
replacing one of the k vertices by a vertical edge with dual 
weight X (denoted with parentheses): 

(x) 



y 



We will also need the following theorem. 

Theorem 4.5. The edge operators Aj^y and Bj^x, using all 
x,y ^ 1, generate PSL(y (k)) (overW) when n>4. 

We will prove Theorem 14.51 in a later paper lfl9ll . In the 
meantime, we include a Sage code that can establish Theo- 
rem |43] for ^ = 8 and for any specific numerical n ifisll . By 
the proof of Proposition l3.2l this is enough to establish that the 
model is universal for linear computation and therefore #P- 
hard to estimate. See also Aharonov, Arad, Eban, and Landau 
for closely related results 

By analogy with Proposition l3.2l Proposition l4.3l and The- 
orem l4.5l together imply that linear circuits can be encoded as 
graphs in the Potts model when n > 4 and y ^ 1, in the same 
way that quantum circuits can be encoded as links using the 
Jones polynomial at a root of unity. 

Thus, Aaronson's theorem again impUes that multiplicative 
approximation of the Potts value Z{G,y,n) is #P-hard when 
« > 4, y < 1, and G is planar This is the first part of Theo- 
remO 

For the second part of Theorem |1.3l it has long been known 
ITtIi that the Jones polynomial V{L,t) of alternating links L 
is equivalent up to normalization to the Potts model Z{G,x,y) 
or the Tutte polynomial with x = —t and y = — f^' (or vice 
versa) for planar graphs G. One argument for this is that by 
replacing each edge of G with a crossing, we can convert a 
contraction-deletion formula for the Potts model to the Kauff- 
man relations. Since then « = f + 2 + f^', we can conclude 
that multiplicative approximation of the Jones polynomial is 
#P-hard when t is real and positive and t ^ I. 

Remark. The Jones polynomial for arbitrary links is equiva- 
lent to a certain multivariate, planar Potts model, and the state 
spaces of the Jones polynomial are isomorphic to those of the 
Potts model. We have already stated that the dimension of 
either state space is a Catalan number 

5. FINAL REMARKS AND QUESTIONS 
5.1. Approximation witli extra information 



In this case the other k—l vertices are upper and lower bound- 
ary, while the edge connects an upper vertex to a lower vertex. 
There are k positions for both of the operators, to make oper- 
ators Aj y and Bj^ for I < j < k; call them edge operators. 

In general, a set of edge operators generates a semigroup 
rather than a group. However, the calculus of compositions 
implies that 

-1 



Therefore: 



Proposition 4.4. Consider the edge operators Ajy and Bj ^ 
generated by a set of weights Y and the associated dual 
weights X. IfX and Y are both closed with respect to taking 
reciprocals, then these edge operators generate a subgroup of 
the projectivized linear group PSL(V'(A;)). 



Theorem ll.2l savs that value-independent approximation of 
certain values of the Jones polynomial are #P-hard even when 
the link L is taken to be a knot. We conjecture that L could 
in addition be a prime knot or even an atoroidal knot. Maybe 
other such restrictions on the structure of L could be imposed. 
But without a result such as that distinguishing the unknot 
(say) is hard, it is not feasible to add arbitrary interesting topo- 
logical restrictions on L to Theorem ll.2l Maybe there is a fast 
algorithm to distinguish the unknot. The Jones polynomial 
would then be easy to compute for knots that are recognized 
as the unknot or as some other specific knot. 



5.2. Other kinds of approximation 

There are many other kinds of partial information about the 
Jones polynomial without any interesting complexity bound 
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to our knowledge. Is the degree of the Jones polynomial in- 
tractable? Is it intractable to determine when the Jones poly- 
nomial vanishes? What if the Jones polynomial is reduced 
mod p for some prime pi 



5.3. Denseness may be more than necessary 

It is easiest to see that a set of gates is universal for lin- 
ear computation if they densely generate an appropriate Lie 
group. For instance, they might generate PSL(2",C) if they 
act on n qubits, or PSL(2",M) inside it. However, dense gen- 
eration is more than necessary for certain types of universal- 
ity. For example, 2-qubit gates that are in PSL(4,Z) do not 
densely generate any Lie group. But suppose that the gate 



/4 -3 
3 4 


1 0\ 
1 


1 

Vo 1 



0/ 



is given |0} for its left qubit and that its left output is also 
postselected to |0). Then the gate reduces 




which is proportional to a rotation by an irrational angle, and 
which does densely generate PSL(4,R) together with other 
integer gates. We do not know the right criteria on linear gates 
to establish #P-hardness results. 



5.4. Morse algorithms may be optimal 

It is common practice to compute the Jones polynomial by 
a strategy known variously as a Morse algorithm, dynamic 
programming, a scanline algorithm, or a divide-and-conquer 



algorithm. (Morse theory in geometric topology is a theory 
of analyzing a topological object by dividing it into horizon- 
tal slices.) For a knot in a plat presentation, the strategy is 
to numerically compute the action of the braid group on the 
skein space. This type of algorithm requires simple exponen- 
tial time and space in the number of strands of the braid, or for 
other kinds of knot diagrams, the width of the diagram. This 
is much better than a direct recursive evaluation of the Jones 
polynomial using a finite set of skein relations; the time com- 
plexity of any such direct algorithm is instead exponential in 
the number of crossings. 

It is natural to wonder whether there are other clever al- 
gorithms that can compute the Jones polynomial even faster. 
However, the proof of Theorem 11.21 could be evidence that 
Morse algorithms are essentially optimal for many kinds of 
knot diagrams. In short, if braids are evaluated using the Jones 
polynomial at the dense roots of unity of Theorem 11.21 then 
they are a model of general planar quantum circuits. 

In more detail, consider a typical hard search problem based 
on classical circuits, and an analogous problem based on 
quantum circuits. For instance, let (z,w) = C{x,y) be a re- 
versible circuit whose input {x,y) and output (z,w) are each 
divided into two registers of equal length. Then it is N P-hard 
to determine whether there is a solution to (z,0) = C{x,Q). 
We conjecture that there are linear-depth, planar circuits C for 
which this problem requires exponential time in |x|. 

Moreover, using denseness at a convenient root of unity and 
the Solovay-Kitaev theorem, this circuit problem can be en- 
coded in a braid with polylogarithmic overhead. We conjec- 
ture that extra logarithmic factors are not necessary for hard- 
ness. We conjecture that there exist linear-length plat presen- 
tations of knots, such that the Jones polynomial requires expo- 
nential time in the bridge number b to estimate at a dense root 
of unity. Such conjectures are very difficult to prove uncondi- 
tionally, because they would imply that #P is not contained in 
FP. But if they were accepted, then Morse algorithms would 
be accepted as essentially optimal. 
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